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SUMMARY 

Most published work relating to dynamically loaded journal bearings are 
directed to determining the minimum film thickness from the predicted journal 
trajectories. These do not give any Information about the subsynchronous whirl 
stability of journal bearing systems since they do not consider the equations 
of motion. It Is, however, necessary to know whether the bearing system opera- 
tion Is stable or not under such an operating condition. 

The purpose of the present paper Is to analyze the stability characteris- 
tics of the system. A linearized perturbation theory about the equilibrium 
point can predict the threshold of stability; however It does not Indicate 
postwhirl orbit detail. The linearized method may Indicate that a bearing Is 
unstable for a given operating condition whereas the nonlinear analysis may 
Indicate that It forms a stable limit cycle. For this reason, a nonlinear 
transient analysis of a rigid rotor supported on oil journal bearings under 
(1) a unidirectional constant load, (2) a unidirectional periodic load, and 
(3) variable rotating load are performed. 

In this paper, the hydrodynamic forces are calculated after solving the 
time-dependent Reynolds equation by a finite difference method with a succes- 
sive overrelaxation scheme. Using these forces, equations of motion are solved 
by the fourth-order Runge-Kutta method to predict the transient behavior of the 
rotor. With the aid of a high-speed digital computer and graphics, the journal 
trajectories are obtained for several different operating conditions. 


INTRODUCTION 

There are two principal approaches by means of which one can analyze the 
whirl instability of a rotor supported on fluid film bearings. These are: 

(1) linearized (perturbation) method and (2) nonlinear transient analysis. In 
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the linearized method a small perturbation of the journal center about Its 
equilibrium point Is given. The stiffness and damping coefficients are deter- 
mined after solving the basic differential equation. These coefficients are 
then used In the equations of motion to find the critical mass parameter and 
whirl ratio. The mass parameter, a function of rotor speed, is a measure of 
stability. The nonlinear transient analysis, on the other hand, gives the 
journal locus and from this one can know about the system stability. 

Stability analysis of finite journal bearings by the linearized method has 
been given by Allaire (ref. 1), whereas Akers, Mlchaelson, and Cameron (ref. 2) 
studied the same bearing configurations using a nonlinear transient approach. 

In reference 2 It was shown that under certain operating conditions the journal 
motion was bounded and could form a limit cycle. 

The aim of the present paper Is to study theoretically the stability char- 
acteristics of finite oil journal bearings under dynamic load using a nonlinear 
transient method. A few papers (refs. 3 to 5) deal with the dynamically loaded 
bearings to predict the journal locus. As these do not consider the equations 
of motion for the prediction of the position of journal center, they cannot 
Indicate whether the bearing system Is stable or not. However, these are use- 
ful for estimating minimum film thickness of dynamically loaded bearings. The 
dynamical equations of motion are solved by fourth-order Runge-Kutta method to 
find eccentricity ratio and attitude angle and their derivatives for the next 
time step. These values are then Introduced In the two-dimensional time- 
dependent Reynolds equation to find the hydrodynamic forces. Following the 
above approach a nonlinear transient analysis of a rigid rotor on oil journal 
bearings under (1) a unidirectional constant load, (2) a unidirectional peri- 
odic load, and (3) variable rotating load Is performed. A number of trajecto- 
ries have been obtained with the aid of a high-speed digital computer and 
graphics . 
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NOMENCLATURE 

radial clearance 
journal diameter 
eccentricity 
hydrodynamic forces 

F r = Fj-C^/nwR^L, Fq = FqC^/owR^L (dimensionless) 

film thickness, Fi = h/C (dimensionless) 

mass parameter, M = MCu 2 /W 0 dimensionless) 

film pressure, p = pC 2 /nc*>R 3 L 

journal radius 

dimensionless time, T = co p t 


2 



t 

x, z, e, z 
H, W 
Wo* Wq 
e 

n 


* 1 , *2 


®1» 

Q 

O) 


Wp 


time 

coordinates, 0 » x/R, z = z/(L/2) (dimensionless) 
load, W => WC 2 / n wR 3 L (dimensionless) 
steady-state load, W 0 = WC 2 /nwR 3 L (dimensionless) 
eccentricity ratio, e » e/C (dimensionless) 
absolute viscosity of oil 
attitude angle 

angular coordinates at which film commences and cavitates, respec- 
tively, from minimum film thickness 

0] ■ ir + , 02 ■ ir + i|>2 

whirl ratio, Q = «p/u> (dimensionless) 

angular velocity of journal 

angular velocity of whirl 


THEORY 


The basic differential equation for pressure distribution In the bearing 
clearance under dynamic conditions can be written as (See fig. 1.) 

h ( h3 I?) + fi ( h3 fi) = 6nR ( w - 2 ft) 15 + 12n H? (1 

Equation £1) when nondimens iona_n zed with the following substitutions: 

0 = x/R, z = z/(L/2), h = h/C, p = pC 2 /n«)R 2 , T = u p t, and Q = up/u, will read 



6 (' - 2Q If) I * ,2Q If 


The film thickness (dimensionless) is 


( 2 ) 


h = 1 + c cos 0 


with the use of equation (3), equation (2) can be written as 


3_ 

90 



6(1 - 2Qcp)(-e sin 0) + 120e COS 0 


(3) 


( 4 ) 


where <p = 9<p/9T and e = 3e/9T. In this study we have used the Reynolds 
boundary conditions, which are given by: 
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p<e,±n - o 


( 0 , 0 ) = 0 


( 5 ) 


and 


p<e,z) = o 


©2 < © < ©i 


91 

30 


(0, z) - 0 


where 0i and ©2 are the angular coordinates at which the film commences and 
reforms, respectively. 


Equation (4) Is solved numerically for pressure by a finite difference 
method with a successive overtaxation scheme satisfying the above boundary 
conditions. Initially, e and <p are set equal to zero to obtain the steady- 
state hydrodynamic forces F r and Fq. 


The forces are computed from 


and 


where 




cos 0 d0 dz 




sin 0 d0 dz 


(6) 


F r C‘ 


nuR^L 


and 


F C 
r 0 u 


6 n uR 3 L 


Releaslng the journal from the steady-state position, one can compute e, 
<p, e, and cp for the next time step by solving the following equations of 
motion: 


MC 


d 2 e 


MC 


dt 

j 2 


2 ‘ \dty 


= F r + W cos <p 


dt 


r* <&)( m ) 


= F q - W sin <p 


The dimensionless form of equations (7) and (8) are: 


Mfi 2 e - flQ 2 e<p 2 - — - COS cp = 0 
W_ 


(7) 

( 8 ) 

(9) 
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(10) 


MQ^e + 2Rfl^e<p - ^-+ Sin <p * 0 


where 


m MCu 

"■if 


and W„ 


W o C 


nwR^L 


The steady-state load W 0 Is obtained by letting e and <p equal to zero. 

Equations (9) and (10) are second-order differential equations In e and 
<p. These are solved by using a fourth-order Runge-Kutta method for constant 
values of £), M, F r , and F©. 

In the following section three types of load are considered. 


A Unidirectional Constant Load 

Assuming the Initial conditions given In table I, the hydrodynamic forces 
under steady-state condition are found. Equations of motion are solved to 
obtain e, e, cp, and <j> for the subsequent time step. Now the new value of 
e, e, and <p are Introduced In equations (3) and (4) to determine the hydro- 
dynamic forces. These forces along with the steady-state load, mass parameter 
and whirl ratio are utilized for the solution of equations (9) and (10). The 
process Is repeated until we get a trajectory that describes the status of the 
system. 


A Unidirectional Periodic Load 
The applied load was assumed to be 

W =. W Q [l + slnQ t)] (11) 

where W 0 Is the applied steady-state load consistent with the eccentricity, 
e 0 = 0.8. At each time steq a new value of W was calculated from equa- 
tion (11). For values of M and 0 given In table I and using W, F r , and 
Fq, equations (9) and (10) were solved for e, e, <p, and <j> for the next time 
step. The rest of the procedure Is repeated as described. 


Variable Rotating Load 

This type of loading Is Interesting since It pertains to an engine bear- 
ing. The data used for the analysis refer to the Ruston and Hornsby 6 VEB-X 
MKIII connecting rod bearing (ref. 5) and Is listed In table I. The polar load 
diagram Is represented by figure 2. For the magnitude of load at different 
positions of crank angle the reader may refer to appendix 4.1 of the paper 
(ref. 5). The time step AT - ir/18 was taken In this case to match with the 
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given applied load at 10° crank angle Interval. The resultant applied load was 
nondimensional Ized using the expression W * WC^/nwR^L. This load was Intro- 
duced In the equations of motion to determine the position of the journal cen- 
ter (e, e, cp, and <p) . The hydrodynamic forces F r and F@ were, however, 
computed after solving equation (4). 


RESULTS AND DISCUSSION 

To assure ourselves that our formulation Is consistent with reference 2, 
our results are shown In figure 3(a) and compared to those of reference 2 shown 
In figure 3(b). The comparison Indicates very good agreement for the condi- 
tions stated in the figure. In figure 4(a) a typical three-dimensional pres- 
sure distribution for a particular time step Is shown. From this figure one 
can see the cavltated region and the variation of pressure throughout the com- 
plete clearance space of the bearing. Figure 5 gives the journal locus when 
the bearing Is stable under the action of constant unidirectional load. In 
figure 3, the journal trajectory within the clearance circle Is shown for an 
unstable bearing. The journal locus reaches a limit cycle. This type of 
observation has been made by many workers (refs. 6 to 8) while dealing with 
short and Infinitely long bearings. 

Having compared the result of the present solution with that of Akers 
et al . (ref. 2), an attempt was made to study the effect of periodic and varia- 
ble load. Figure 6 shows the journal locus for a periodic load superimposed on 
the constant load given the same as those of a unidirectional constant load. 
From these two figures It may be seen that a bearing which Is stable under the 
action of a unidirectional constant load can be made unstable when a periodic 
load is applied. In the latter case the journal locus reaches a limit cycle. 

The journal center trajectory of the connecting rod bearing (variable 
rotating load) is shown in figure 7. The trajectory Is very complex: it nei- 

ther tends to go to an equilibrium point nor to reach a limit cycle. It may be 
mentioned that the value of the minimum film thickness for this case exceeds 
that obtained by others using the mobility method of solution. The mobility 
method (ref. 3), however, does not consider the stability of the system. 

The present solution does not take account of the transport of fluid 
through the cavltated region and consequently does not take full account of the 
oil film history. A subsequent analysis has been undertaken that will Include 
the oil film history effects for these dynamic conditions and will be compared 
at a later date. A proper accounting of the mass flow In time Is believed to 
be Important in dynamically loaded bearings and should Influence the trajectory 
of the rotor in a more realistic way. 


CONCLUSIONS 

From the nonlinear transient analysis of an oil-film journal bearing under 
different dynamic loads with Reynolds type boundary conditions, the following 
conclusions are evident. 
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1. Although the analysis Is costly In terms of computer time, It gives the 
orbital trajectory within the clearance circle which Is not obtainable using a 
simplified (linearized) theory. 

2. As It Is shown by others while dealing with short and Infinitely long 
bearing theory, the journal locus for a finite bearing ends In a limit cycle 
for an unstable system. 

3. A stable journal with a constant unidirection load can be made unstable 
by superimposing a periodic load on the system. 

4. A more correct (cavltated) boundary condition using the prehistory of 
film may show some Interesting phenomena which are not revealed In the present 
method of solution. 


REFERENCES 

1. Allaire, P.E., Design of Journal Bearings for High-Speed Rotating 
Machinery, Fundamentals of the Design of Fluid Film Bearings, ASME 
Publication, New York, 1979, pp. 45-83. 

2. Akers, A., Mlchaelson, S., and Cameron, A., Stability Contours for a 
Whirling Finite Journal Bearing, J. Lub. Tech., Trans. ASME, Series F, 
vol . 93, no. 1, 1971, pp. 177-190. 

3. Booker, J.F., Dynamically Loaded Journal Bearings - Mobility Method of 
Solution, J. Basic Engr. Trans. ASME, Series D, vol. 187, no. 3, 1965, 
p. 537. 

4. Horsnell, R. and McCall Ion, H., Prediction of Some Journal Bearing 
Characteristics Linder Static and Dynamic Loading, Proc. Lubrication and 
Wear Convention, Inst. Mech. Engrs., London, 1963, pp. 126-138. 

5. Campbell, J. Love, P.P., Martin, F.A., and Reflque, S.C., Bearings for 

Reciprocating Machinery: A review of the Present State of Theoretical 

Experiment and Service Knowledge, Conf. Lubrication and Wear, Proc. Inst. 
Mech. Engrs., vol. 182, (Part 3A), 1967-68, pp. 51-74. 

6. Kirk, R. and Guenter, E.J., Transient Journal Bearing Analysis, NASA 
CR-1549, June 1970. 

7. Holmes, R., The Vibration of a Rigid Shaft on Short Sleeve Bearings, J. 
Mech. Engr. Sc 1 . , vol. 2, no. 4, 1960, p. 337. 

8. Jakobsen, K. and Christensen, H., Nonlinear Transient Vibrations In 
Journal Bearings, Proc. Inst. Mech. Engr. vol. 183, Part 3P, 1968-69, 
pp. 50-56. 


7 


TABLE I. - INPUT DATA FOR VARIOUS CONDITIONS OF LOAD 


Condi ti on 

L/D 

Unidi rectional 
constant 

1.0 

Unidi rectional 
periodi c 

1.0 

Variable 

rotating 

0.562 


M 

Q 

5. 

0.5 

5. 

0.5 

5. 

1 .0 
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RELATIVE TO CYLINDER AXIS 
FIGURE 2. - POLAR LOAD DIAGRAM OF ENGINE BEARING. 



(B) OBTAINED BY AKERS ET AL. 123 


FIGURE 3. - JOURNAL CENTER TRAJECTORY. (L/D = 1.0, 
M = 10.24, 0 = 0.5). 






FIGURE <1. - A TYPICAL PRESSURE DISTRIBUTION (L/D = 1.0, £ 0 ■ 0.8, 
M = 5, 0 = 0.5, T - 0). 



FIGURE 6. - JOURNAL CENTER TRAJECTORY FOR A UNIDIRECTIONAL PERIODIC 
LOAD (L/D = 1.0, £ n = 0.8, M = 5, 0 = 0.5). 


10 






NASA 

National Aeronautics and 
Space Administration 


Report Documentation Page 


1. Report No. NASA X M- 102309 

AVSCOM TR 87-C-26 


2. Government Accession No. 


3. Recipient’s Catalog No. 


4. Title and Subtitle 

Stability of a Rigid Rotor Supported on Oil-Film Journal 
Bearings Under Dynamic Load 


5. Report Date 


6. Performing Organization Code 


7. Author(s) 

B.C Majumdar and D.E. Brewe 


8. Performing Organization Report No. 
E-3727 


9. Performing Organization Name and Address 

NASA Lewis Research Center 
Cleveland, Ohio 44135-3191 
and 

Propulsion Directorate 

U.S. Army Aviation Research and Technology Activity— AVSCOM 
Cleveland, Ohio 44135-3127 


10. Work Unit No. 

505-63-81 


11. Contract or Grant No. 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 

Washington, D.C. 20546-0001 

and 

U.S. Army Aviation Systems Command 
St. Louis, Mo. 63120-1798 


13. Type of Report and Period Covered 
Technical Memorandum 


14. Sponsoring Agency Code 


15. Supplementary Notes 

Prepared for the National Seminar on Bearings Madras, India, September 17-18, 1987. B.C. Majumdar, National Research Council— 
NASA Research Associate; present address: Indian Institute of Technology, Department of Mechanical Engineering, Kharagpur-721302, 
India. D.E. Brewe, Propulsion Directorate, U.S. Army Aviation Research and Technology Activity— AVSCOM. 


16. Abstract 

Most published work relating to dynamically loaded journal bearings are directed to determining the minimum film 
thickness from the predicted journal trajectories. These do not give any information about the subsynchronous 
whirl stability of journal bearing systems since they do not consider the equations of motion. It is, however, 
necessary to know whether the bearing system operation is stable or not under such an operating condition. The 
purpose of the present paper is to analyze the stability characteristics of the system. A linearized perturbation 
theory about the equilibrium point can predict the threshold of stability; however it does not indicate postwhirl 
orbit detail. The linearized method may indicate that a bearing is unstable for a given operating condition whereas 
the nonlinear analysis may indicate that it forms a stable limit cycle. For this reason, a nonlinear transient 
analysis of a rigid rotor supported on oil journal bearings under (1) a unidirectional constant load, (2) a unidirec- 
tional periodic load, and (3) variable rotating load are performed. In x this paper, the hydrodynamic forces are 
calculated after solving the time-dependent Reynolds equation by a finite difference method with a successive 
overrelaxation scheme. Using these forces, equations of motion are solved by the fourth-order Runge-Kutta 
method to predict the transient behavior of the rotor. With the aid of a high-speed digital computer and graphics, 
the journal trajectories are obtained for several different operating conditions. 


17. Key Words (Suggested by Author(s)) 

Journal bearings; Stability; Dynamic loads; 
Bearings; Dynamics 


18. Distribution Statement 

Unclassified - Unlimited 
Subject Category 34 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21 . No of pages 

22. Price* 

Unclassified 

Unclassified 

12 

1. AQ3 1 


NASA FORM 1626 OCT 86 


*For sale by the National Technical Information Service, Springfield, Virginia 22161 




